Effect of the accelerating growth of communications networks on their structure 
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Motivated by data on the evolution of the Internet and World Wide Web we consider scenarios 
of self-organization of the nonlinearly growing networks into free-scale structures. We find that 
the accelerating growth of the networks establishes their structure. For the growing networks with 
preferential linking and increasing density of links, two scenarios are possible. In one of them, the 
value of the exponent 7 of the connectivity distribution is between 3/2 and 2. In the other, 7 > 2 
and the distribution is necessarily non-stationary. 
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The incredible place of the Internet in our civilization 
evokes the exponentially increasing flow of the studies 
of evolving networks (Internet, World Wide Web, neural 
networks, nets of citations, collaboration networks, etc.) 
Mm ■ The study of communications networks has a long 
history (e.g., see the papers of "the father" of the Inter- 
net, P. Baran jlQ]). Nevertheless, the first data on the 
structure of the Internet, Web, etc were obtained only 
recently B 0,O ]l3[ . Most of them are on the simplest 
"one-site" characteristic of networks - the connectivity 
distribution. The observation of power-law dependences 
of the connectivity distribution in these networks puts 
forward the problem: how a growing network can orga- 
nize itself into a scale-free structure? 

The most natural mechanism of such self-organization 
is the network growth with the preferential attachment 
of new links to sites with a high number of connections 
P, [lll[l4| . There are several ways to introduce the prefer- 
ential linking which lead to different values of scaling ex- 
ponents of the_ connectivity distribution inside of a huge 
interval ll 



Nevertheless, it is still unclear why the 
exponents of the real systems have their specific values, 
and why the connectivity distributions have their specific 
forms? E.g., the last measured value for the scaling ex- 
ponent of the distribution of a number of incoming links 
in the World Wide Web is 2.1 [|jl|Jl|. 

The last data on the Internet evolution demonstrate 
that the total number of links increases more quickly than 
the number of sites ||l2[. For instance, the study of the 
inter-domain topology of the Internet ll^ has shown that 
the relation of the edges and the nodes equaled 1.71 in 
November of 1997 (3015 nodes) and 1.88 in December of 
1998 (8256 nodes). The number of pages in the World 
Wide Web was 203 x 10*5 with 1466 x 10^ links between 
them in May 1999 (the corresponding ratio is 7.22) and 
it was already 271 x 10^ pages and 2130 x 10^ finks in 
October 1999 (the ratio is 7.86) ^. That looks quite 
naturally - old sites may establish new connections all 
the time. Old documents of the Web may be updated 
and new references may be added to them. Collabora- 
tion networks becomes to be more dense while growing 
because of increasing possibility to find a collaborator. 



Thus we face certainly the wide spreaded situation: the 
density of the links of a network becomes to be more and 
more high during its evolution, the total number of links 
grows more fast than the number of sites, the growth is 
nonlinear, that is accelerating. 

This very significant factor is omitted in most of con- 
sidered models p, pT| , p^ -17|. One may ask: What effect 
has the accelerating growth on the structure of the net- 
works? What kind of the acceleration can produce scale- 
free networks? What form can take the connectivity dis- 
tribution in such networks? 

In the present Communication we consider these prob- 
lems using general arguments. Assuming that the net- 
works are scale-free, we describe the possible connectiv- 
ity distributions and show that, in this case, the total 
number of links is a power-law function of the network 
size. We introduce simple models with preferential link- 
ing accounting for the accelerated growth and study their 
evolution both analytically using the continuous connec- 
tivity approach ||lj,|l^,|l^ and by simulation. The arising 
distributions cover the range of the forms predicted from 
the general consideration. Also, we find the situation in 
which the fiow of new links, increasing via a power law, 
takes out the growing network from the class of free-scale 
structures. 

We start from the most general considerations. In the 
scale-free networks, a wide range of the connectivity dis- 
tribution function is of the power-law form, P{q) oc q~'^ . 
It will be clear from the following that, to keep the net- 
work in the class of free-scale nets, the fiow of new links 
has to be a power function of the number of sites of the 
network that plays the role of time, i.e., be proportional 
toi". 

First, let us assume that the exponent of the distribu- 
tion is less than two. The reasonable range is 1 < 7 < 2. 
To produce the restricted average connectivity (that is 
proportional to i"), the distribution has to have a cut- 
off at large q. Its natural value is of the order of the 
total number of sites, that is proportional to t""*"^. It is 
also clear from the following that, if the distribution is 
non-stationary, the connectivity distribution of the scale- 



free networks has to be of the form, t' 



q 



between 



q r^ t^ and q ~ t°'~^^, where z and x are some exponents. 
It is a restricted function in the range below q ^ t^. 
z > since the network grows. From the normahza- 
tion condition, J^ dqP{q,t) = 1, we get x — zjia — !)■ 
The average connectivity, g is of the order t"+^/t, then 
r - f*' dqqt^'q-'f - ^z+{2-i)ia+i) f^^^ie value of 
the integral is defined by its upper limit). Therefore, 
a = z + {2--f){a + l), and, finally, 7 = l + (l + 2:)/(l + a). 
One sees that 2 < a to keep the exponent 7 below 2 as 
it was assumed. Also, one sees that the lower bound- 
ary for 7, 1 + 1/(1 + a), is approached for the stationary 
distribution, z ^ 0. In this case, the form of the distribu- 
tion is completely fixed by the accelerating growth, the 
exponent 7 depends only on a. 

We do not know the values of a for real networks. Ob- 
viously, a cannot exceed 1 (the total number of links has 
to be smaller than i^/2 since one may forbid multiple 
links). Hence, 7 > 3/2. The value a — corresponds 
to permanent density of links ||l^. The density of con- 
nections in the real networks remains rather low all the 
time 112] , so one may reasonably assume that a is small. 
Therefore, the lower boundary of the possible values of 7 
is close to 2. 

The other possibility is 7 > 2. In this case, again 
X — 2/(7 — 1) but the integral for the average connectiv- 
ity is defined by its lower limit, t" ~ J^^(-,-i) dqqV'q^'^ ^ 
t^-^T-2)/(T-i). Hence, 7 = 1 -h z/a, and z > a. Thus, 
we have described the possible forms of the connectivity 
distribution. 

Let us demonstrate how these distributions may arise 
in the nonlinearly growing networks with the preferen- 
tial linking. We do not restrict ourselves to some model 
generating the increasing fiow of new links dynamically 
in the process of the growth since it would narrow the 
range of possible dependencies on time for this flow. (We 
do not discuss here possible mechanisms of such gener- 
ation. One of them was considered in |18|.) Instead of 
that, we prefer to introduce this dependence ab initio. 

Let us introduce the following simple model of the 
growing network with directed links. We consider only 
the distribution of incoming links, so, here, the connec- 
tivity is the number of incoming links qs of the site s. Let 
each instant, a new site be added to the network. Let it 
initial connectivity be n, where n > 0. That means that 
it has n incoming links from some non-specified old sites. 
Let extra coi" new directed links be distributed between 
old sites. We assume that each of these links comes out 
from an unspecified site and is directed to some site s 
chosen with probability proportional to the sum of its 
connectivity and some constant, qs + A. Here n + A has 
to be positive. 

In fact, we generalize the model |1^ to the case of the 
increasing number of new links. The introduced power- 
law dependence of the number of new links keeps the 
network in the class of free-scale networks but changes 
crucially the connectivity distribution. 

Note that, unlike the first model of a scale- free growing 



network with preferential linking proposed by Barabasi 
and Albert [^, the constant A is introduced in the rule 
of the distribution of the new links between the sites of 
the network. The sum n + A characterizes the initial 
"attractiveness" of a site for new links |13| . 

We use the continuous connectivity approach |14| , [l5| , [l7| 
which gives exact values for the scaling exponents of such 
systems as it was demonstrated in |15|] . Each site is la- 
beled by the time of its birth, < s < i. The connectivity 
distribution of the site s in the continuous-g approxima- 
tion has the ioTm.p{q, s,t) = S{q — q(s,t)) fla, where S{ ) 
is the delta-function and q{s, t) is the average connectiv- 
ity of the site s at time t. 

The resulting equation for q{s,t) is very simple: 



dq{s,t) 
dt 



cot" 



qis,t)+A 



jQdu[q{u,t)+A] 



(1) 



with the boundary condition q(t, t) — n. One may also 
set g(0,0) — 0. The equation describes the distribution 
of flow of new links between sites according the intro- 
duced rules. Applying /^ ds to both sides of Eq. (|l|) one 
gets /g duq{u,t) — nt + cot"^^/{a + 1), so 



cot" 



dt {n + A)t + coi«+V(a + 1) 

The general solution of Eq. (Ill) is 



[qi.s,t)+A]. (2) 



qis,t)+A^g{s) 



1 
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1 



l+a 



1+1/c 



Co 



t 



-(a-l/a) 



(3) 



where g{s) is an arbitrary function of s. It is fixed by the 
boundary condition, q(t, t) = n. Thus one obtains the 
result for g(s, t), 



q{s,t)+A 
n + A 



l + {n + A)il + a)t-°'/co 
l + {n + A){l + a)s-'^/co 



i-i-i/( 



(i) 



In the interval [{n + A){1 + a)/co]^/" < s < t, 

-ia+l) 



q{s,t)^{n + A) (0 



-(«+i) 



(4) 



(5) 



Thus, the exponent /3, q{s,t) oc s ^ , equals l + a and is 
bigger than 1. The dependence q{s) becomes constant, 



g(s,f) = (n + A)-VW-££- 



1+1/0 



+Q+1 



(6) 



at s <C [(n + A){1 + a)/coY^°'. One may compare the 
result, Eq. (0), with the total number of links in the 
network, N{t) « coi"+7(l + ")• 

In the continuous-g approximation, one may easily find 
the distribution P{q, t) using the obtained q[s, t): 



P{q,t) = j [ dsp{q,S,t) = ^ [ 
^ JO ^ Jo 



ds6{q — q{s,t)) 



1 /dq{s,t)\-' 



t \ ds 



(7) 



< 



Therefore, in the region 1 <C q/(ri + A) 
{ co/[{n + A)il + a)] ji+^Z" fi+"^ the connectivity distri- 
bution has the foUowing form: 

^ (n + A)V(i+") ^_[i+i/(i+„)] 
1 + a 

Thus, we obtain the stationary connectivity distribution, 
P{q) oc q~^ with 7 = 1 + 1/(1 + a) that befongs to one 
of the types described above. 

The last result follows also from the relation /3(7— 1) = 
1 between the scaling exponents, which may be obtained 
from the assumption that q{s,t) and P{q, t) are power- 
law functions of s and q correspondingly p5[ . Note that 
HA -^ 00, the network is out of the class of scale- free 
networks for any a > 0. 

The introduction of the increasing flow of new links 
in the problem changes crucially the distribution of the 
connectivities. Indeed, in the case of a constant density 
of links, the value of the scaling exponent 7 of the con- 
nectivity distribution varies from 2 to infinity depending 
on the network parameters, n and A p5| , and the values 
of the exponent /? are between and 1. Here, for the 
increasing density of links, we obtain 7 below 2 and the 
exponent (3 exceeds 1. The values of the exponents are 
independent of n and A. 

To demonstrate the other possibility, 7 > 2, we con- 
sider below the model with a different rule of the distribu- 
tion of new links. We make the only change in the studied 
above model. Let now a new link be directed to some site 
s with probability proportional to qs{t)/q{t) + B, where 
qs{t) is the connectivity of site s, q{t) is an average con- 
nectivity of sites of the network, and B is some posi- 
tive constant. q{t) = coi"/(l -I- a), so the new links are 
distributed between sites with probability proportional 
to qs(t) + BcQt"/{l + a), where ag is positive constant. 
Hence, the previously introduced A becomes to be time- 
dependent. 

Note that, if a = 0, this rule for the preferential link- 
ing is also reduced to the known model of the growing 
network with a permanent input flow of new links |15[| . 

Repeating the previous calculations, one gets the equa- 
tion 



dq{s,t) 
dt 



cof" 



q{s,t)+Bcot°'/{l + a) 



nt + Bcof^+^/il + a) + Cf)t°'+^ / {a + 1) 

(9) 



with the boundary condition q{t,t) — n. At long times, 
one obtains 



dq{s, t) l + a q{s, t) + Bcq i"/(l + a) 



dt 



l + B 



t 



(10) 



The solution of Eq. (hO) is 



q{s,t) = 



Bcqs" 



1- Ba 



SX -(l+a)/(l+S) Bcot° 



(l) 



1 -Ba 



(11) 



If B = 0, we obtain the previous result, /? = l + a. For 
s">n(l-Ba)/(Bco), 



q{s,t) 



Boot'' 



l-Ba \\t 



a-(l+a)/{l+B) 



1 



(12) 



Therefore, the scaling exponents of the growing network 
are (3 = {1 + a)/{l + B) - a = (1 - Ba)/{\ + B) 
and 7 = 1 -h 1//3 - 1 + [(1 + a)/(l + B) - a]-^ = 
2 + B{1 + a)/(l — Ba). The connectivity distribution 
differs sharply from the distribution obtained for the 
previous model. It is nonstationary and is of the form 

P{q,t) ^ i-l + (l+")(l-Sa)g-[l+(l+B)/(l-Ba)] fo^ ^ ^ ^a^ 

In this case, (3 < 1 and 7 > 2 for any positive a and 
B. The scaling regime is realized if Ba < 1. The general 
phase diagram for both considered model is shown in Fig. 

Note that, in both considered cases, one cannot set 
a = directly in the obtained expression for the scaling 
exponents. In such a situation, we get from Eqs. ( p[) o r 
(|) /3 = [1 + (A + 7i)/co]-i and 7 = 2 + (A + n)/col|. 

It is known that the used continuous approach gives 
exact results for the scaling exponents of the growing net- 
works with a constant density of connections |15| . Nev- 
ertheless, it is approximate, so we have checked the ob- 
tained above results by simulation. 

The results of the simulation of considered models are 
shown in Figs. ^ and y. The size of networks in both 
studied cases is 10000 sites. The number of the at- 
tempts equals 1000. In Fig. g, we present the log-log 
plots of the average connectivity vs number of a site for 
a — 0.5, 71 = 1, A = 1.0, Co = 1.0 (the flrst model) and for 
a = 0.5, n = 1,5 = 0.15, Co = 1.0 (the second one). In 
Fig. ^, for these values of parameters of the models, we 
show the log-log plots of the connectivity distribution. 

The obtained values of the scaling exponents are within 
the error of the simulation from the corresponding ones 
found analytically. The values (3 = 1.46 (1.5) are ob- 
tained from the simulation and analytically (in brack- 
ets) for the first model with the written out parameters, 
(3 — 0.85 (0.804) are the corresponding values for the sec- 
ond model. 7 = 1.69 (1.667) and 7 = 2.19 (2.243) are the 
values of the critical exponent of the connectivity distri- 
bution obtained for the first and for the second models, 
relatively. One may see that the correspondence is really 
good. 

Several different values of the scaling exponent of the 
distribution of incoming links in the World Wide Web 
were published (as far as we know, any data on the ex- 
ponent j3 are absent yet). The available data are 7 = 2.1 
§,|l]jlD and 1.94 |l9). The strong difference between 



these values may be explained by the difference of the 
sizes of the scanned areas of Web. The most huge area 
was studied in |1^], so the value 7 = 2.1 seems to be the 
best one. As we have noted, one may assume reasonably 
that a is small in the real networks. We have shown that, 
in such a situation, the lower boundary for the possible 
values of 7 is slightly below 2. We have demonstrated 
that, for 7 > 2, the connectivity distribution has to be 
non-stationary if the growth of the network is acceler- 
ating. There are no data that let us learn whether the 
connectivity distributions of the World Wide Web and 
the Internet are stationary or not. Our results make this 
question intriguing. 

The World Wide Web is still in the initial stage of its 
evolution. Perhaps, the parameters of the accelerating 
growth will change. In this case, our answers demon- 
strate the possibility of changing of 7. We have shown 
that it may become even less than 2 in future. 

To demonstrate all the existing possibilities we have 
considered the models of growing networks with the par- 
ticular rules of the preferential attachment of new links. 
These models cover the range of possibilities but provide 
us only with particular values of the scaling exponents. 
Of course, there exists a lot of additional factors (aging 
pTf and dying |0,|l^ of sites, etc.) which may change 
these particular values. 

In summary, we have studied the nonlinear, accelerat- 
ing growth of the scale- free networks. We have demon- 
strated that it can be one of the most significant factors 
defining their structure. We have described the possi- 
ble connectivity distributions of such networks and have 
fixed the lower boundary for the scaling exponent 7. 
Only the power-law time-dependence of the input flow 
of new links can keep the network inside of the class 
of scale-free networks. Nevertheless, we have found the 
region of parameters in which the scale-free structure is 
impossible (see the phase diagram in Fig. |l|). Our re- 
sults demonstrate possibility of quite different scenarios 
of the network evolution and let us hope to approach 
satisfactory description of the real networks. 
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FIG. 1. Phase diagram of the networks with the accel- 
erating growth under consideration. The networks are out 
of the class of scale-free nets ( "exponential" ) above the line 
a = 1/B. The exponent 7 equals 3 on the dashed line and 
2.1 (value for the World Wide Web) on the dash-dotted one. 
7 < 2 on the line 5 = 0. 



FIG. 2. Log-log plot of the average connectivity of a site FIG. 3. Log-log plot of the distribution of the number 
vs its number (birth time) for the considered models. For the of incoming links of sites for the considered models. For the 
first model, a — 0.5, n — 1,A — 1.0, co — 1.0. For the second first model, a — 0.5, n = 1,A — 1.0, co = 1.0. For the sec- 
model, a — 0.5, n = 1,B = 0.15, cq — 1.0. The dashed lines ond model, a = 0.5, n = 1,_B = 0.15, cq = 1.0. The dashed 
have the slopes equal to the values of the scaling exponent /3 lines have the slopes equal to the values of the scaling ex- 
obtained analytically. ponent 7 obtained analytically. For better presentation, the 

dependences are displaced along the vertical axis. 
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